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ABSTRACT 

The 10/zm silicate feature is an essential diagnostic of dust-grain growth and 
planet formation in young circumstellar disks. The Spitzer Space Telescope has 
revolutionized the study of this feature, but due to its small (85cm) aperture, 
it cannot spatially resolve small/medium separation binaries (<3"; < 420 AU) 
at the distances of the nearest star-forming regions (~140 pc). Large, 6-10m 
ground-based telescopes with mid-infrared instruments can resolve these systems. 

In this paper, we spatially resolve the 0.88" binary, UY Aur, with 
MMTAO/BLINC-MIRAC4 mid-infrared spectroscopy. We then compare our 
spectra to Spitzer/IRS (unresolved) spectroscopy, and resolved images from 
IRTF/MIRAC2, Keck/OSCIR and Gemini/Michelle, which were taken over the 
past decade. We find that UY Aur A has extremely pristine, ISM-like grains 
and that UY Aur B has an unusually shaped silicate feature, which is probably 
the result of blended emission and absorption from foreground extinction in its 
disk. We also find evidence for variability in both UY Aur A and UY Aur B by 
comparing synthetic photometry from our spectra with resolved imaging from 
previous epochs. The photometric variability of UY Aur A could be an indica- 
tion that the silicate emission itself is variable, as was recently found in EX Lupi. 
Otherwise, the thermal continuum is variable, and either the ISM-like dust has 
never evolved, or it is being replenished, perhaps by UY Aur's circumbinary disk. 
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Introduction 



Grain growth, which is one of the initial steps to planet formation, can be studied by 
observing the 10/zm silicate features of young stars. By modeling the shape and strength 
of the feature, it is often possible to determine the different dust grain constituents (amor- 



(Bouwman et al. 


2001 


Sareent et al. 


2009|). ) 



2009 



metrics to determine the characteristic grain-size and crystallinity of dust (van Boekel et al. 



g 

20031 : iPrzvgodda et al.ll2003l : lApai et alEoOd iKessler-Silacci et al.ll2006l : iPascucci et al.ll2008h . 
Based on this morphological classification, a logical evolutionary sequence involv es small, 



amor phous, ISM-like grains combining and annealing into large (or more porous — see 
20081 ) and crystalline g rains on the way to forming planetesimals and eventually planets. For 
a complete review, see iNatta et al.l (120071 ) . 



Voshchinnikov fc Hem 



The Spitzer Space Teles cope and its I RS sp ectrograph have observed numerous young 
stellar objects, including the iFurlan et al.l (120061 ) Taurus- Aurigae survey of 111 Class II and 
III objects. However, a large fraction of t hese results ar e conta minated by unresolved binaries 
due to Spitzer's small (85cm) aperture. IFurlan et al.l (120061 ) have noted that along with a 
large number of "normal" silicate spectra, they have 15 outlier spectra, some of which may 
be unresolved blends. Known binaries exist in both the normal subsample and the outlier 
subsample. Complementing the space- based results is a spa tially-resolved, N-band (~8/mi- 
13.5/xm) imaging survey of binaries by iMcCabe et al.l (120061 ) . However, there is still a need 
for resolved spectroscopy of these binaries. 

In this paper, we present resolved, broad N-band spectroscopy of the 0.88" binary UY 
Aur using the MMT deformable secondary adaptive optics system (MMTAO), and the mid- 
infrared camera, BLINC-MIRAC4. We also present unpublished resolved imaging data from 



1 The observations reported here were partially obtained at the MMT Observatory, a facility operated 
jointly by the Smithsonian Institution and the University of Arizona. 

2 The observations reported here were partially obtained at the Gemini Observatory, which is operated 
by the Association of Universities for Research in Astronomy, Inc., under a cooperative agreement with the 
NSF on behalf of the Gemini partnership: the National Science Foundation (United States), the Science 
and Technology Facilities Council (United Kingdom), the National Research Council (Canada), CONICYT 
(Chile), the Australian Research Council (Australia), CNPq (Brazil), and CONICET (Argentina) 

3 The observations reported here were partially obtained at the Infrared Telescope Facility, which is 
operated by the University of Hawaii under Cooperative Agreement no. NCC 5-538 with the National 
Aeronautics and Space Administration, Science Mission Directorate, Planetary Astronomy Program. 
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two earlier epochs using IRTF/MIRA C2 and Gemini /Mic helle and compare our results to 
published Keck/OSCIR pho t ometry (iMcCabe et al.l 120061) and Spitze r/IRS blended spec- 



troscopy (IFurlan et al.l 120061 ; IWatson et al 



200 



9; 



Pascucci et al.l 120081 ) . Collectively, these 



observations span more than a decade. 

UY Aur is a particularly interesting T Tauri system in the Tauru s- Aurigae star forming 
region that was originally discovered to be a binary in the optical (jJoy fc van Biesbroeck 
19441 ). before one of its compo nents, UY Aur B, faded to more than 5 magnitudes fainter 
than the primary at R-band Jrlerbst et al.l \l 9m. The sy s tem w as observed to be a bi- 
nary in the infrared by iGhez et al. I Jl993h and iLeinert et al.l (119931 ) , and subsequent H-band 
and K-band imaging reveal e d that the system is highly variable at those wavelengths too 



flKenvon fc Hartmannl Il995l : IClose et al.l Il998l : IWhite fc Ghed 1200 ll ; iBrandeker et all 12003 



Hioki et al. l2007h . The source of the variability, while assumed to be m ostly variable ac- 



cretio n/extinction around UY Aur B may come from UY Aur A a s well (IBrandeker et al. 



200 31). However, s patially resolved, absolute calibration data from IClose et al.l (119981 ) and 
Hioki et al.l ( 120071 ) show that UY Aur A's H-band flux stayed nearly constant (AH =0.4 
magnitudes) for three data points between 1996 and 2005, while UY Aur B varied by a 
maximum of AH =1.3 magnitudes. 



Dutrey et al.l ( 119961 ) and iDuvert et al.l ( 119981 ) resolved a Keplerian circumbinary disk 
around UY Aur A+B usin g 13 CO millimeter interferometry. The disk was subsequently 



resolv ed in scattered light by lClose et al.l (119981 ) and pol arized scattered lig ht by lPotter et al. 
( 120001 ) . Using scattered-light multi- wavelength data, IClose et al.l (119981 ) showed that the 



dust composition of the circumbinary disk was dominated by small (0.03/zm-0.6/im) grains 
and postulated, based on their images, that spiral-like streamers from the circumbinary 
disk were replenishing small grains in the circumstellar disks of the binary components. 
The circumstellar disks are thought to be truncated to <l/2 the se paration of the binary 
(<60 AU assuming a face-on orbit) based on dynamical simulations (lArtymowicz fc Lubow 
1994). The system also has a bipolar outflow, although which star is its source is unknown 
jHirth et al.lll997h . 



Hartigan &: Kenyonl (120031 ) used the Space Telescope Imaging Spectrograph (STIS) to 



acquire optical spectra of both binary components. Their results show that UY Aur A is 
an M0 classical T Tauri star with a mass of 0.6 M© and a fairly high accretion rate of 
10~ 7,64 M & /yr. UY Aur B was characterized as an M2.5 classical T Tauri star with a mass 
of 0.34 M & and a similarly high accretion rate of 10~ 7 ' 70 M Q /yr Extinction measurements 
showed that UY Aur A was extincted by ^4^=0.55 and UY Aur B was extincted by Ay=2.65. 
This large disparity in the extinction of Class II binaries is uncommon since both stars lie 
along a common line of sight in the dark cloud. Possible causes of this phenomenon are 
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extinction of UY Aur B from a foreground object (such as the disk of UY Aur A), or self- 
extinction from the disk (caused, for example, by UY Aur B's disk potentially being viewed 
from close to edge-on). 

In summary, UY Aur is a binary classical T Tauri system with two truncated circum- 
stellar disks. UY Aur B is an infrared companion that has exhibited variability in the optical 
and near-IR. UY Aur A may be variable in the near-IR as well. The system is surrounded 
by a circumbinary disk, which may contain spiral-like streamers to the circumstellar disks. 



Of the 111 young stars in lFurlan et al.l (120061 ) 's sample, UY Aur an d LkCa 15 haye spec - 
tra that are the best examples of pristine, ISM-like silicate features (IWatson et al.l 120091 ). 
which could be the result of a lack of dust-grain evolution or variable and/or replenishing 
mechanisms. Given that these Spitzer measurements are of a blended spectru m, and UY Aur 
B is h eavily extincted at optical and near-IR wavelengths, it is surprising that IWatson et al. 
(120091 ) would measure such a pristine (de-extincted) ISM-like feature. Our resolved spec- 
troscopy will help determine the nature and differences of the dust components in UY Aur 
A and B. 



2. Observations and Reductions 



2.1. MMTAO /BLINC-MIRAC4 AO Spectra 

On January 14, 2009 UT, we observed UY Aur with MMTAO /BLINC-MIRAC4 using 



Llovd-Hart 


2000; 


Wildi et al.ll2003: 


Brusa et al. 



20041 ). BLINC-MIRAC4 (Bracewell Infrared Nulling Cryostat; Mid-IR Array Camera, Gen. 
4) is a combined AO-optimized mid-IR imager (MIRAC4) and nulling interferometer (BLINC), 
which for these observations, is used in its "imaging" mode. MI RAC4 is functionally equiv- 



alent to previous versions of MIRAC (e.g. 

umented in Hinz et al (in prep). BLINC is descibed in iHinz et al 



Hoffmann et~allll998h. with its new features doc- 

(hoooh. U sing AO, 



BLINC-MIRAC4 can achieve Strehls of up to ~98% at 10//m (IClose et al.ll2003af ). 



MIRAC4 was first used in its sp ectroscopic mode for grism trace-centroiding, without 
a flux calibration (ISkemer et al.l 120091 ) . This paper presents the first flux-calibrated spectra 
with MIRAC4, and thus a full description of our observing strategies and reductions will 
follow. 



MIRAC4 operates at a variety of platescales ranging from 0.055"/pixel to 0.110"/pixel. 
In order to observe the full N-band range (7.6/im-13.3/mi), we operate close to the coarsest 
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platescale (0.105"/pixel). For all the spectra in this paper, we use an oversized (1") slit, to 
negate object alignment effects. In this configuration, MIRAC4 has a spectroscopic range of 
7.6/zm-13.3/zm, which is limited by an N-band filter to block higher order grism modes. Our 
wavelength calibration was done in a lab the same week as the observations and resulted in a 
measured dispersion of 27 nm/pixel. No telluric wavelength calibration is done at this time 
because of the lack of available lines in the N-band atmospheric window. Because we use 
an oversized slit, our spectral resolution is limited by diffraction. At the center of N-band 
(10.5/xm), the diffraction limit of the 6.5 meter MMT is ~ X/D =0.32". This corresponds 
to a spectral resolution of ~125. 

Conditions on January 14, 2009 UT were photometric at the beginning of the night 
and we observed the mid-IR standard, e Tau, at two different airmasses in broad-N spectra, 
followed by the science target, UY Aur. Ideally, the UY Aur observations would have been 
followed by another observation of e Tau. However, the wind shifted and strengthened 
after our UY Aur observation, precluding a third observation of the standard. All of our 
observations used adaptive optics with a loop-speed of 550 Hz. Details of our observations 
are presented in Table [TJ 

Each object was aligned in the slit for two nod beams separated by 5" (we chop per- 
pendicular to the slit with an 8" throw). Because we were observing the binary UY Aur, we 
rotated the camera with respect to the sky to manually align the binary with the slit. The 
telescope rotator was turned on to maintain the position angle. 

All three gro ups of spectroscopic observations were reduced with our custom artifact 



removal software ( iSkemer et al.l 1200a ) and mean combined. We compared the two reduced 
e Tau spectra to determine the quality of our flux calibration. Because the airmass of e 
Tau changed between our observations (see Table [T|, we do a telluric airmass correction by 
dividing our spectra by atmospheric transmission curves for Mauna Kea at the appropriate 
airmasses (we assumed a 3 mm water vapor column density but using different values did not 
significantly change our r esults; the t ransmission curves are courtesy of Gemin jj and use the 
ATRAN code described in 



(Ord 



19921 ). This technique has also been described in lHonda et al. 
(l2004j ) and lMcCabe et al.l (120061 ) among others. After telluric airmass correction, the ratio of 
the divided e Tau spectra is shown in Figured], boxcar-smoothed with a 0.19/zm kernal. The 
large feature centered at ~9.5/im is telluric ozone, and it is clear that within the grey-boxed 
region, our flux calibration is poor. We measure the mean and standard deviation of the 
residuals and find that outside the ozone feature, our flux calibration has an error of 2.3% ± 



4 http://www. gemini.edu/sciops/telescopes-and-sites/observing-condition-constraints/transmission- 
spectra 
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3.0%. The mean difference between the spectra (2.3%) is treated as a global error, and the 
standard deviation (3.0%), while not completely uncorrelated with wavelength, is treated as 
a local random error. The same calculation yields an error of 19.6% ± 9.0% error inside 
ozone. This could be caused by an incorrect assumption of ozone density in our telluric 
correction (ozone varies with latitude and season) or a density fluctuation of ozone during 
our observations (fast day/night or temperature/pressure induced fluctuations). The best 
way to address both of these problems is to group future observations closely in both airmass 
and time. 

O ur UY Aur observations were reduced similarly to those of e Tau. At a separation of 
0.88" Jdose et alill998l : lMcCabe et al.ll20o3 l. the binary is easil y resolved by MMTA O /BLINC- 
MIRAC4, which has superresolved binaries as tight as 0.12" (jSkemer et al.ll2008l ). We de- 
termine the relative flux of the two components of UY Aur at every wavelength by fitting 
each 1-D slice with a corresponding 1-D slice from the second (closest in airmass and time) 
observation of the point source e Tau. The images of the spectra of e Tau and UY Aur as 
well as a best-fit residual plot are shown in Figure [2j 



We perform an absolute ca 



ibration on UY A ur A and B by comparing them to the 



G9.5III mid-IR standard, e Tau (jCohen et al.lll999l ) and using the telluric airmass correction 



technique described above in this section. The error analysis on our calibrated spectra is 
somewhat complex due to the different sources and correlations between errors. We list the 
five major error sources here: 



measurement error due to photon noise, slight PSF mismatch, and crosstalk artifacts on 
the MIRAC4 chip. The crosstalk artifacts manifest in our flux measurements as occa- 
sional bad measurements at a single wavelength (detector row). To robustly measure 
the flux of our spectra in a way that removes the effects of an occasional bad single-row 
measurement, we median combine our data in 7-pixel (0.19/im) bins (which effectively 
degrades our spectral resolution by a factor of two). We empirically calculate the dis- 
persion in our median measurements by taking the separation between the second and 
sixth ordered values in each 7-pixel bin, and multiplying by the Monte Carlo simula- 
tion determined value, 0.303, which gives a 1-sigma (gaussian) error bar on the median. 
Typical errors range from l%-4% of the flux with errors up to ~5% of the flux inside 
telluric ozone. 



globa l mo del error (i.e. correlated throughout the full spectrum) on the ICohen et al. 
(119991 ) model of e Tau is 2.259% throughout the spectrum. 



local model error on the ICohen et al 
across the spectral range. 



f|1999h model of e Tau ranges from 1.032%-1.265% 
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global telluric error based on our comparison of the two e Tau measurements shown in 
Figure [TJ We assume a global error of 2.3% outside the ozone feature and 10% inside 
ozone. The 10% error inside ozone is less than the 19.6% discrepancy shown in Figure 
[TJ However, we believe 10% to be more appropriate based on the fact that the second 
e Tau measurement and UY Aur measurement from Table [1] are separated by less in 
airmass/time than the two e Tau measurements. 

local telluric error based on our comparison of the two e Tau measurements shown in 
Figure [TJ We assume errors of 3.0% outside the ozone feature and 5% inside ozone. 
Again, the 5% error inside ozone is less than the 9.0% discrepancy shown in Figure [TJ 
but is reasonable considering the proximity of our observations. 



All five error sources are combined in quadrature and plotted about the median binned 
spectra in Figure |3j The figure shows four spectra: UY Aur A, UY Aur B, the combined UY 
Aur A and B, and a reference Spitzer/IRS spectrum of the unresolved UY Aur binary (Jeroen 
Bouwman and Ilaria Pascucci, private communication). Measurement errors (and total 
errors) are calculated separately for each spectrum. The Spitzer/ IRS spectrum uses Spitzer 



pipeline S18.7.0 and methods described in iBouwman et al.l (120081) . as wel l as a correction for 



pointing offest errors using the methods described in ISwain et al.l (120081 ) . The Spitzer/IRS 
spectrum has been binned to the same wavelength range as the MMTAO/BLINC-MIRAC4 
data, and measurement errors are calculated similarly to the MMTAO /BLINC-MIRAC4 data 
to include the "error" resulting from binning a non-neglible spectral range. The dominant 
source of error for Spitzer is a 3% global model uncertainty. 

The Spitzer/IRS spectrum has slightly more flux than the MMTAO /BLINC-MIRAC4 
UY Aur combined spectrum. However, the difference is not significant enough to infer 
variability based on the overall flux levels. In Section 4, we will show that the shape of the 
spectra and the relative fluxes of the binary components suggest temporal variability. It is 
also evident from this plot that UY Aur A and UY Aur B have similar overall fluxes but very 
different silicate profiles. This will be discussed in Section 4.1 when we analyze the shapes 
of the silicate features. The flux ratio of the two components is plotted as a function of 
wavelength in Figure HJ Because relative photometry has fewer error sources than absolute 
photometry, flux ratios are powerful measurements for determining differences between the 
silicate feature (Section 4.1) and variability (Section 4.2). 
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2.2. Gemini/Michelle Imaging 

On September 27, 2004 UT, Gemini/Michelle observed UY Aur using Michelle's spec- 
troscopic mode for Gemini program GN-2004B-Q59. We obtained acquisition images of UY 
Aur (which is easily resolved) and the standard star t Aur in the Si-5 filter (11.15/im-12.25/im 
at half- maximum). A summary of these observations is presented in Table [TJ We did not 
attempt to use the spectroscopic data as the binary was aligned perpendicular to Michelle's 
slit. The data were reduced with the standard Gemini IRAF pipeline and are shown in 
Figure [51 



We performed photometry on UY Aur with IRAF daophot (jStetsonl 119871 ) using t Aur 
as a PSF a nd photometric standard, i Aur is a K3II Cohen mid-IR standard ([Cohen et al. 
19991 . boOlh . which has been flux-calibrated to 63.074 Jy in the Si-5 filtefl The 



error on our 

absolute flux calibration is dominated by variable sky conditions, since 1 Aur and UY Aur 
were observed 15 minutes apart. Based on the well- matched airmass and weather conditions, 
and the fact that the data were taken on a relatively dry night in a clear region of the N-band 
window, we assign a 5% absolute photometric error to our measurements. The uncertainty 
in our UY Aur A-B flux ratio measurement is most likely dominated by PSF fitting error. 
Between the two measurments of UY Aur, the flux ratio varied by ~1%. We conservatively 
adopt 2% as the error on flux ratio. The results of our Gemini/Michelle imaging photometry 
are presented in Table [2J 



2.3. IRTF/MIRAC2 Imaging 



On November 15, 1998 UT, IRTF/MIRAC2 observed UY Aur using MIRAC2's imaging 
mode with the older (now replaced) 128x128 detector. We obtained images of UY Aur (which 
is marginally resolved) and the standard star a Tau in the 10.3/im filter (9.73/im-10.73//m 
at half-maximum). A summary of these ob servations is present ed in Table [TJ The data were 
reduced with the custom IRAF pipeline of I Close et al.l (l2003bl ) and are shown in Figure [6j 



We performed PSF fitting photometry on UY Aur with IRAF daophot ([Stetson! 119871 ) 
using a Tau as a PSF and photometric standard, a Tau is a MIRAC2 standard, w hich 
has been flux-calibrated to 584 Jy in the 10.3 /xm filtefl using the I Cohen et al.l (I1992al ) sys- 
tem defined by a Lyr. The night was photometric and both the object and the standard 
were observed at very low airmass (1.02-1.09). However, the 10.3/im filter partially over- 



J http://www. gemini.edu/sciops/instruments/mir/MIRStdFluxes.txt 
6 cfa- www. harvard.edu/~jhora/mirac/mrcman.pdf 
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laps with the telluric ozone feature (which can adversely effect absolute calibrations) so we 
conservatively assume an absolute calibration error of 8%, which is added in quadrature to 
the photon/fitting errors. We also adopt 4% as the error on the flux ratio from the errors 
given by the IRAF allstar task. The results of our IRTF/MIRAC2 imaging photometry are 
presented in Table [2j 



3.1. 



3. Analysis 
Normalization of the Spectra 



Following the example of numerous authors ( 


van Boekel et al. 


2003; 


Przveodda et al. 


2003; 


Kessler-Silacci et al. 


2006; 


Pascucci et al. 


2008), we continuum subtract and normalize 



our spectra (see Figure [7]) in order to measure size/shape parameters of the 10/xm silicate 
feature, which can be used to infer the characteristic dust grain size, and crystallinity. After 
continuum subtraction, the normalization follows the formula 



S v = l + 



F - F 

v - 1 u,c 

< F v ,. > 



where F v is the spectrum and F V)C is the fit continuum and the denominator, < F U)C > is the 
continuum averaged over the wavelength range of the silicate feature. 

We perform a Monte Carlo simulation of the error terms described in Section 2.1 to 
construct 10000 sets of synthetic spectra. The error terms are correlated where appropriate 
so that errors from the absolute calibration (global and local model errors and global and 
local telluric errors) are the same for the three MMTAO/BLINC-MIRAC4 spectra (UY Aur 
A, UY Aur B, and combined UY Aur) and global model/telluric/calibration uncertainties 
are the same throughout each full spectrum in a given Monte Carlo trial. For each trial 
we perform a best-fit linear continuum subtraction based on the continuum at 8.0/xm and 
12.5/iHi-13.0/im. Since the N-band atmospheric window precludes us from measuring the 
quadratic and higher order continuum terms, we adopt the quadratic term of the continuum 
from the Spitzer/IRS spectrum (fit between 5.3/mi-8.0/zm and 12.6/im-14.2 /xm; to avoid 



silicat e) for all of our spectra. As can be seen from an independent analysis by lWatson et al. 



(120091 ). the higher order (non-linear) continuum terms for UY Aur are small. The Monte 
Carlo analysis gives us a range of values for the norm alized flux at each wavelength from 
which we calculate error bars. iPascucci et al.l (120081 ) use a similar Monte Carlo analysis 
(with different error sources). The 1-a range of the fit continua is shown in Figure [7] along 
with the normalized fluxes derived from equation [1] and our Monte Carlo analysis. A few 
of the normalized spectra are also plotted against each other for ease of comparison in 
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Figure EJ Table [3] shows the 9.9/iin normalized flux and the ratio of the 11.25/xm/9.9/zm 
normalized fluxes, which are used as a proxy for grain-size and crystallinity respectively. We 
also measure the ratio of these parameters for UY Aur A and B (since this helps supress much 
of the uncertainty caused by the absolute error terms). We find that the ratio of the 9.9/im 
normalized flux terms is 1.40±0.06 and the ratio of the 11.25/zm/9.9/zm normalized color is 
0.82±0.04. Because the quadratic term of our continuum subtraction is not well known (we 
adopt it from a different epoch and it is not resolved for the two binary components), we 
may have some residual systematic errors that we do not present. However, the difference 
between using the quadratic term from the Spitzer data and using no quadratic term at all 
creates a discrepancy of only 0.01 for the values in Table [3] so we believe this effect to be 
small. 



We compare the photometry of multiple epochs of UY Aur observations to check 
for photometric variability, presenting both absolute photometry of the full system and 
relative photometry of the resolved binary for the ground-based data in Table [2j With 
the spectroscopic data, we can construct synthetic filters to compare our actual photom- 
etry (IRTF/MIRAC2, Keck/OSCIR, Gemini/Michelle) to the spectra (MMTAO /BLINC- 
MIRAC4 and Spitzer/IRS). For example, photometry in the Gemini/Michelle Si-5 filter is 
compared to MMTAO/BLINC-MIRAC4 spectra that have been convolved onto the Gem- 
ini/Michelle Si-5 photometric system. This necessitates a thorough understanding and cali- 
bration of each photometric system. 

The traditionally quoted energy flux of an object through a photometric system is given 



where F\(\) is the spectrum of the object as a function of A and -R(A) is the system response 
function (i.e. the product of the sky transmission, the system QE and the filter profile as a 
function of A). If we convert this to an F v system (i.e. Janskys), the equation becomes 



3.2. Photometric Variability 



by 



F(filter) 



J F\(\) ■ R(X)d\ 
jR{X)d\ 



(2) 



F {filter) 



f F V (X) ■ c ■ A" 2 ■ R{X)dX 
jR{X)dX 



(3) 



However, when comparing detector counts f or relative and abso lute calibration, it is necessary 
to integrate photons rather than energy (jBessell et al.lll998l ). which necessitates that we 
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divide the flux term in equation [3] by E = hv. This gives us 

Count Sl _ J F U<1 (X) ■ \- 1 -R(X)d\ 
Counts 2 ~ S F Vi2 {X) ■ A- 1 ■ R(\)d\ 



(4) 



where countsi is the ratio of photons detected for object 1 compared to object 2, with fluxes 

countS2 1 j ? 

F Ut i and F Vj2 respectively. 

In this formalism, Equations |2] or [3] should be used to determine absolute fluxes for 
photometry (i.e. to put a Cohen model spectrum onto a photometric system) while Equation 
H] should be used for photometry with actual data (i.e. counts measured for object 1 versus 
counts measured for object 2 in a photometric system). 

In order to put the Spitzer/IRS and MMTAO /BLINC-MIRAC4 spectra onto the pho- 
tometr ic systems of IRTF/ M IRAC2, Ke ck/OSCIR and Gemini/Michelle, we use Equation 
Hwith lCohen et all jl992afi . Il995l Il999h spectra of a Tau (for IRTF/MIRAC2), a Lyr (for 
Keck/OSCIR) and i Aur (for Gemini/Michelle) as F v2 {X). We use the telluric absorption 
spectrum at 1.0 airmass es for Mau na Kea assuming 1.0 mm precipitable water vapor as 



described in Section 2.1 (|Lordlll992h . and we use filter curves for 10.3/int| (IRTF/MIRAC2), 
N-bancfl (Keck/OSCIR) and Si-£@ (Gemini/Michelle). 

Equation [Hand its corresponding uncertainty is calculated using a Monte Carlo approach 
with the same error sources (on F u i and F v ^ 2 ) and functionality as the analysis described in 
Section 3.1. The binary flux ratios and corresponding uncertainties are calculated using the 
relative errors shown in Figure HI 



4. Discussion 

4.1. Normalized Spectra 

Our continuum-subtracted, normalized spectra in Figure [7] show that UY Aur A's 
silicate emission mostly dominates the combined normalized spectrum, but that UY Aur 
B's role is not insignificant. A direct comparison of the MMTAO/BLINC-MIRAC4 nor- 
malized spectrum (Figure [S]) with the Spitzer normalized spectrum shows that the fea- 
ture itself has changed, as the Spitzer spectrum has a stronger feature that could be in- 



7 http : / / www. cfa. harvard . edu/ ~ emama j ek / mirac / index, html 
8 http://www. gemini.edu/sciops/instruments/oscir/oscirlndex. html 
9 http://www. gemini.edu/sciops/instruments/michelle/imaging/filters 
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terpreted as smaller, less-evolved grains. Interestingly, the bottom panel shows that the 
MMTAO/BLINC-MIRAC4 normalized spectrum of UY Aur A is consistent with the Spitzer 
normalized spectrum of the unresolved pair. This could mean that UY Aur A's feature was 
even more dominant at the time of the Spitzer observations. Given that we have observed 
significant variability in both binary components (see below in Section 4.2), it is possible 
that changes in either star's silicate feature could have caused this. Without two epochs of 
resolved spectra (the Spitzer spectrum is a blend), it is impossible to know if UY Aur A's 
silicate feature has varied, or just its continuum. 



As was noted by I Watson et al.l (120091 ) . the Spitzer spectrum of the unresolved UY Aur 
binary shows almost uniquely clear evidence for small, pristine, ISM-like grains. In Figure |9l 
we compare our normalized spectrum of UY Aur A with th e scale d optical depth of the ISM 
towards Sgr A* using the measurements of iKemper et al.l (120041 ) who find that the ISM is 
dominated by ~0.1/xm amorphous grains. Our spectrum is consistent with these ISM grains 
with a reduced xl value of 0.75. Given that UY Aur A is variable (see below in section 4.2) 
and has pristine ISM-like grains, it will be interesting to see if future resolv ed spectroscopy 



finds a varying silicate feature similar to the recent observations of EX Lupi ( 1 Abraham et al. 
2009h . 



While it would be interesting to compare the dust properties of UY Aur A and B, 
extinction in the spectrum of B precludes a direct comparison. Instead, we can de-extinct 
the spectrum of UY Aur B, to see if we can reproduce UY Aur A's ISM-like silicate feature 
with a plausible amount of extinction. We use the ratio of the fluxes between UY Aur A and 
B (Figure HJ which has lower errors than our flux-calibr a ted sp ectra) and de-extinct UY Aur 
B using scaled ISM optical depths from IKemper et al.l (|2004|). These scaled optical depths 
can be converted to Ay using the ISM extinction curve of iRieke fc Lebofskyl ( 119851 ). With 
this, we do a 3 parameter best-fit (Ay and a linear continuum) to match the spectrum of 
UY Aur B to UY Aur A. Our results are plotted in Figure [Till where the black diamonds are 
the observed flux ratios between UY Aur A and B (Figure H]), and the teal asterisks are the 
theoretical flux ratios between UY Aur A and B after a best-fit de-extinction of UY Aur B 
plus the addition of a linear continuum. The best-fit extinction value is Ay = 5.1 with a 
reduced \ 2 V of 1.71. This is a plausible amount of extinction given that UY Aur B has been 
observed t o vary by >4.5 mag nitudes at R-band, due to a mixture of variable extinction and 
accretion (IHerbst et al.lll995l ). Figure [TT1 shows the normalized spectru m of UY Aur B after 
being de-extincted by Ay = (as observed), 2.65 (the value measured by lHartigan fc Kenyon 
20031 ) and 5.1 (the best-fit value to make UY Aur B's silicate feature like UY Aur A's). If 
extinction in UY Aur B is variable, then 5.1 magnitudes of Ay extinction would allow the 
dust grains in UY Aur A and B to be similar. If extinction in UY Aur B is not variable 
or was lower than Ay=5.1 at the time of our observations, as was true when measured 
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by lHartigan fc Kenyonl (120031 ). then the silicate feature in UY Aur B is shallower than A's 
which would mean it has larger grains. Disentangling these two scenarios is only possible 
with simultaneous resolved extinction measurements and resolved mid-IR spectroscopy. 



4.2. Photometric Variability 



As is evident from Table [2] and a corresponding plot shown in Figure [121 the UY Aur 
system is variable in the mid-infrared. Among our results, we note that UY Aur A and 
UY Aur B became significantly fainter in the Si-5 filter between September 27, 2004 and 
January 14, 2009. Also, between November 16, 1999 and January 14, 2009, UY Aur A got 
fainter in N-band while UY Aur B got brighter. This change is further highlighted in the 
flux ratio measurements which dropped from 2.07±0.03 to 1.19±0.01 over that period. Our 
results are somewhat surprising considerin g that UY Aur B was a ssumed to be the dominant 
source of variability in the near-infrared (IBrandeker et al.l 120031 ) . The fact that UY Aur A 
is also variable in the mid-infrared may mean that it deserves a more complete study at 
shorter wavelengths as well. Currently, the best resolve d near-infrared var iabili ty data set. 



with a bsolute calibrations, is 3 H-band measurements by Idose et al.l (Il998h and iHioki et al. 
( 120071 ) that show ~0.4 magnitudes of variability for UY Aur A and ~1.3 magnitudes of 



variability for UY Aur B between 1996 and 2005. If this variability is caused entirely by 
extinction, we can use the iRieke fc Lebofskyi ( 119851 ) ISM extinction law to say that N-band 
extinction variability would be ~0.1 magnitudes for UY Aur A, and ~0.4 magnitudes for 
UY Aur B. Based on our mid-IR variability data from Table [21 this hypothetical extinction 
variability would be enough to explain UY Aur B's mid-IR behavior, but not UY Aur A's. 
Thus, some other process, such as variable accretion, must be responsible for UY Aur A's 
mid-IR variability. 



If I Close et al.l (jl998h 's hypothesis (following lArtymowicz &: Lubowl Il996l ) that spiral 
streamers are responsible for dust-replenishment and variability in the system, then our 
mid-infrared data suggest the streamers may be affecting the primary star (UY Aur A) as 
well as the secondary. Future studies should address whether binaries are more likely than 
single stars to be variable in the mid-infrared, and whether the variability and the unevolved 
ISM-like grains are contingent on large circumbinary disks feeding truncated circumstellar 
disks, such as the ones surrounding UY Aur and GG Tau. Spectroscopic measurements, 
which are able to determine if the source of the variability is continuum or optically thin 
dust emission, would be particularly useful. 
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5. Conclusions 

While the Spitzer IRS spectrograph has observed numerous young stellar objects to 
study grain-growth, Spitzer's 85cm aperture cannot resolve most binaries (tighter than ~3"). 
As a result, large aperture (ground-based mid-IR systems and/or JWST) must be used to 
separate tight/medium-separation binaries. 

Using the 6.5 meter MMT with MMTAO /BLINC-MIRAC4, we obtained high-Strehl, 
spatially resolved, low-resolution, N-band spectra of the 0.88" binary UY Aur. We compare 
our results to resolved photometry from IRTF/MIRAC2, Gemini/Michelle, and Keck/OSCIR 
as well as an unresolved spectrum from Spitzer. We find the following: 

1) UY Aur A and UY Aur B have very different silicate features, which when combined, 
produce a blended spectrum that includes emission from both sources and absorption from 
UY Aur B. Other binaries may have similarly complicated silicate blends and should be 
resolved by large-aperture telescopes. 

2) UY Aur A's silicate spectrum is consistent with dust emission from pristine, ISM-like 
grains. UY Aur B, which is an infrared companion, has a much flatter silicate spectrum that 
may be the result of foreground extinction from an edge-on disk. As a result, we cannot 
determine whether UY Aur A and UY Aur B share similar dust properties (although if 
corrected for Ay =5.1 magnitudes of extinction, UY Aur B has a similar silicate feature to 
UY Aur A). The fact that UY Aur A exhibits ISM-like dust implies that either its dust has 
not evolved (via grain-growth) or that its small grains are replenished by variable processes 
(such as streamers from the circumbinary disk reservoir). 

3) While UY Aur B has shown optical and near-IR variability in the past, we show that 
UY Aur A and UY Aur B both show significant variability in the mid-IR. The variability 
in UY Aur B is expected based on its presumably changing extinction profile. However, 
the mid-IR variability of UY Aur A, especially considering its unevolved, ISM-like grains, is 
more surprising. Future resolved spectra will be able to determine if this variability is due 
to continuu m changes, or if it is d ue to episodes of dust-processing as was recently observed 



in EX Lupi (Abraham et all 120091 ) . 
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Table 1. Imaging and Spectroscopic Observations of UY Aur 



Date (UT) 


Telescope /Instrument 


Object 


Filter 


Exposure Length (s) 


# of Exposures 


Airmass 


Jan 14, 2009 


MMTAO /BLINC-MIRAC4 a 


e Tau 


Grism 


10 


10 


1.37-1.38 


Jan 14, 2009 


MMTAO /BLINC-MIRAC4 a 


e Tau 


Grism 


10 


14 


1.15-1.16 


Jan 14, 2009 


MMTAO /BLINC-MIRAC4 a 


UY Aur 


Grism 


10 


18 


1.09-1.11 


Sep 27, 2004 


Gemini/Michcllc b 


i Aur 


Si-5 


4 


1 


1.05 


Sep 27, 2004 


Gemini/Michelle b 


UY Aur 


Si-5 


4 


2 


1.03 


Nov 15, 1998 


IRTF/MIRAC2 


a Tau 


10.3Mm 


5 


1 


1.09 


Nov 15, 1998 


IRTF/MIRAC2 


UY Aur 


10.3^m 


10 


54 


1.02-1.05 



a All data were taken with MMTAO at a loop-speed of 550 Hz 

b The data were all taken within a period of 15 minutes and had IQ=70%, CC=50% and WV=20% (Gemini observing 
constraints for Mauna Kea — www.gemini.edu/sciops/telescopes-and-sites/observing-condition-constraints) 
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Table 2. Photometry of UY Aur 



Date 


Telescope / 


Filter 


UY Aur 


UY Aur A 


UY Aur B 


Binary flux ratio 




Instrument 




Flux (Jy) 


Flux (Jy) 


Flux (Jy) 


(A/B) 


Nov 15, 1998 


IRTF/MIRAC2 


10.3/jm 


2.54±0.24 


1.54±0.17 


1.00±0.14 


1.55±0.06 


Feb 27, 2004 


Spitzer/IRS 


10.3/im 


2.81±0.08 








Jan 14, 2009 


MMTAO /BLINC-MIRAC4 


10.3/im 


2.63±0.10 


1.49±0.06 


1.14±0.04 


1.30±0.01 


Nov 16, 1999 


Keck/OSCIR 


N-band 


2.27±0.22 


1.53±0.14 


0.74±0.05 


2.07±0.03 


Feb 27, 2004 


Spitzer/IRS 


N-band 


2.35±0.07 








Jan 14, 2009 


MMTAO /BLINC-MIRAC4 


N-band 


2.17±0.08 


1.17±0.05 


0.99±0.04 


1.19±0.01 


Feb 27, 2004 


Spitzer/IRS 


Si-5 


2.62±0.08 








Sep 27, 2004 


Gemini /Michelle 


Si-5 


3.26±0.16 


1.95±0.10 


1.31±0.07 


1.49±0.03 


Jan 14, 2009 


MMTAO/BLINC-MIRAC4 


Si-5 


2.42±0.09 


1.32±0.05 


1.10±0.04 


1.22±0.02 



Note. — The Spitzer/IRS and MMTAO/BLINC-MIRAC4 data are spectra that have been converted to the photometric 
systems of 10.3/im (IRTF/MIRAC2), N-band ( Keck/OSCIR) a nd Si- 5 (Gemini/Michelle) as described in Section 3.2. The 
Keck/OSCIR photometry of UY Aur is from M cCab e et al Their data is presented in magnitudes, which we 

convert to Jy using the zero-mag flux of 33.99 Jy for a Lyr presented in their paper. 
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Table 3. Silicate Shape/Strength Measurements 



Spectrum 


59.9 


Sn. 25/59.9 


UY Aur (Spitzcr unresolved) 


1.53±0.01 


0.86±0.01 


UY Aur (MMTAO combined) 


1.36±0.05 


0.83±0.04 


UY Aur A (MMTAO) 


1.59±0.07 


0.75±0.04 


UY Aur B (MMTAO) 


1.14±0.06 


0.92±0.06 
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8 9 10 11 12 

Wavelength (/im) 



Fig. 1. — The divided, processed spectra of two observations of e Tau demonstrates the 
quality of our absolute calibration. After a telluric airmass correction (described in text), 
we are able to achieve a calibration discrepancy of 2.3% ± 3.0% outside of the grey-boxed 
region coincident with telluric ozone. Inside ozone the discrepancy is 19.6% ±9.0%. 
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Best-Fit Residuals 

i 



Fig. 2. — Spectra of our PSF/flux calibrator, e Tau, our science target, UY Aur, and the 
best-fit residuals after fitting for the flux of each component of UY Aur. Note: For clarity, 
the e Tau spectrum is presented at 1/10 the stretch of the UY Aur spectra and residuals. 
UY Aur A is 0.88" from UY Aur B. 
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Fig. 3.— MMTAO /BLINC-MIRAC4 flux calibrated spectra of UY Aur A, UY Aur B, their 
combined spectrum and a comparison Spitzer (unresolved) spectrum. The error bars shown 
are the combination of several error sources (described in text). 
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Fig. A.— The flux ratio of UY Aur as a function of wavelength from the MMTAO/BLINC- 
MIRAC4 spectra. Because the relative errors on each components flux are significantly 
lower than the absolute flux calibration errors, flux ratios provide a powerful measurement 
for demonstrating variability and differences between the silicate features. 



- 23 - 



lota Aur (PSR 


UY Aur B UY Aur A 

\ / 





# • 



Fig. 5. — Gemini/Michelle acquisition images of Iota Aur (the PSF/flux calibrator) and UY 
Aur. For Iota Aur, North is up and East is left. For UY Aur, the coordinate system has 
been rotated clockwise by 132 degrees. 
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Fig. 6.— IRTF/MIRAC2 images of UY Aur and Alpha Tau (the PSF/flux-calibrator) along 
with the best-fit residuals of the PSF-fitting photometry. The best-fit residuals are at the 
photon noise floor of our image, which gives confidence that the binary was accurately 
resolved and extracted. For all three images, North is up and East is left. 
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Fig. 7.— LEFT: MMTAO /BLINC-MIRAC4 spectra shown in Figure [3] with the l-a range 
of our Monte Carlo continuum plotted in grey. RIGHT: Normalized spectra (see Equation 
[I]) after continuum subtraction. The error bars are the result of our Monte Carlo analysis 
(described in Section 3.1). 
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Fig. 8. — Comparisons of normalized spectra from Figure [71 The negligible error bars on the 
Spitzer normalized spectra (shown in Figure [7J are removed for clarity. Our combined UY 
Aur normalized spectrum shows statistically significant variability compared to the Spitzer 
spectrum. However, UY Aur A itself is statistically indistinguishable from the Spitzer spec- 
trum. 




Fig. 9. — The continuum-subtracted, normalized spectrum of UY Aur A is compared to the 
(scaled) optical depth of the interstellar medium towards Sgr A*. The curves are consistent 
with reduced \ 2 = 0.75. The fact that UY Aur A's dust is very similar to the ISM either 
means that it has not evolved (via grain-growth), or that variable processes (such as accretion 
from the circumbinary disk) are replenishing the small ~0.1/zm grains. 
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Fig. 10. — In order to see if UY Aur B's silicate feature is similar to UY Aur A's but 
with more extinction and a different underlying continuum, we take the observed relative 
fluxes between UY Aur A and B and do a best-fit de-extinction of B, plus the addition of a 
linear continuum, to try to produce an equal flux-ratio across all wavelengths. The observed 
relative fluxes between UY Aur A and B are shown as black diamonds (as in Figure H]). The 
best-fit de-extincted ratio (+continuum) is shown as teal asterisks and has a reduced xl °f 
1.71 with respect to the dashed line (unity). 
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Fig. 11. — Normalized spectra of UY Aur B after de-extinction compared to the scaled ISM- 
optical depth shown in Figure M Fr om top to bottom, UY Au r B has been de-extincted by 
Ay of 0, 2.65 (the value measured by lHartigan &: Kenyonll2003l ) and 5.1 (the best-fit value to 
make UY Aur B's silicate feature like UY Aur A's). Assuming a large amount of extinction 
towards UY Aur B, it is possible that UY Aur B's silicate grains are similar to UY Aur A's. 
For lower values of extinction, UY Aur B's silicate grains would have to be larger than UY 
Aur A's grains. 
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Fig. 12. — Variability of UY Aur and its two components in three filters. The photometry is 
listed in Table |2j Variability is seen, to some extent, in the combined UY Aur plot, as well 
as in the UY Aur A and B plots. However, ground-based resolved data are most sensitive 
to changes in the flux ratio (seen in the bottom right). Our most recent MMTAO/BL1NC- 
MIRAC4 data show that the binary is much closer to equal flux than it has been in the past 
in all three filters. 
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